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Abstract.  Thermally  excited  plasma  modes  are  observed  in  trapped,  near-thermal-equilibrium  pure 
electron  plasmas  over  a temperature  range  of  0.05  <T  <5  eV.  The  measured  thermal  emission  spec- 
tra together  with  a separate  measurement  of  the  wave  absorption  coefficient  uniquely  determines 
the  temperature.  Alternately,  kinetic  theory  including  the  antenna  geometry  and  the  measured  mode 
damping  (i.e.  spectral  width)  gives  the  plasma  impedance,  obviating  the  reflection  measurement. 
This  non-destructive  temperature  diagnostic  agrees  well  with  standard  diagnostics,  and  may  be  use- 
ful for  expensive  species  such  as  anti-matter. 


Even  in  a stable  plasma  there  are  a finite  level  of  fluctuating  electric  field  and  plasma 
waves.  The  plasma  waves  are  “normal  modes”  [ 1]  of  the  system,  and  represent  degrees  of 
freedom  which  are  excited  in  thermal  equilibrium.  Plasma  waves  are  emitted  by  particles 
as  they  move  about  in  the  plasma,  and  they  are  absorbed  (e.g.  Landau-damped)  by  the 
plasma.  The  balance  between  emission  and  absorption  leads  to  a thermal  level  of  field 
fluctuation.  Thermally  excited  modes  can  be  a diagnostic  tool  for  non-neutral  plasmas  as 
shown  by  Gould  for  cyclotron  waves  [2].  Similarly,  center  of  mass  modes  are  commonly 
used  in  the  “single  particle  regime”  with  a highly  tuned  resonant  circuit  [3]. 

We  describe  measurements  of  spontaneous  thermal  excitation  of  Trivelpiece-Gould 
(TG)  modes.  The  thermal  radiation  observed  comes  from  collective  plasma  modes.  The 
spectrum  of  these  modes  shows  several  well-separated  peaks  at  the  frequency  of  TG 
modes.  The  radiation  we  describe  comes  from  one  of  these  peaks,  typically  form  the 
lowest  longitudinal  and  radial  mode  number  of  an  azimuthally  symmetric  mode. 

Here,  the  spectrum  of  thermally-excited  T-G  modes  is  measured  in  near-thermal- 
equilibrium  pure  electron  plasmas  over  a temperature  range  of  0.05  < kT  < 5 eV,  us- 
ing a room-temperature  receiver.  The  received  spectrum  for  each  mode  is  nominally  a 
Lorentzian  at  frequency  /o  with  width  y,  superimposed  on  a background  noise  level;  non- 
Lorentzian  spectra  also  occur,  when  the  receiver  noise  and  impedance  is  non-negligible. 
By  Nyquist’s  theorem,  these  fluctuations  represent  thermal  noise  generated  by  the  real 
part  of  the  plasma  impedance  Zp  [4].  The  area  of  the  received  spectral  peak  is  propor- 
tional to  kT,  but  it  also  depends  on  the  antenna  coupling  as  represented  by  Zp.  This 
impedance  can  be  measured  directly  by  a separate  wave  reflection/absorption  measure- 
ment, or  it  can  be  calculated  from  a kinetic  theory  of  random  test  particles.  Overall, 
the  technique  allows  a rapid  non-destructive  diagnostic  of  the  plasma  temperature  with 
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FIGURE  1.  Schematic  diagram  of  the  cylindrical  Penning-Malmberg  trap  with  (a)  the  wave  transmis- 
sion/reception electronics;  (b)  the  thermal  emission  equivalent  circuit;  and  (c)  the  reflection/absorption 
electronics. 


±25%  accuracy. 

Fluctuation  measurements  were  obtained  from  pure  electron  plasmas  combined  in  two 
similar  Penning-Malmberg  traps,  “EV”  and  “IV”  [5],  differing  mainly  in  plasma  diam- 
eter and  magnetic  field  strength.  The  EV  trap  consists  of  a series  of  hollow  conducting 
cylinders  of  radius  r„,  = 3.8  cm  contained  in  ultra-high  vacuum  at  P w 10“*'’  Torr  with 
a uniform  axial  magnetic  field  of  B — 375  G.  Electrons  are  injected  from  a hot  tung- 
sten filament,  and  contained  axially  by  voltages  Vc  ^ -200  V on  end  electrodes.  Typical 
plasmas  have  N m 10^  electrons  in  a column  length  Lp  w 24  cm,  with  a plasma  radius 
rp  w 1 .7  cm  and  a central  density  no  « 10^  cm“^.  For  IV,  the  parameters  are  B = 30  kG, 
rp  = 0.2  cm,  r„,  = 2.86  cm,  and  Lp  = 4l  cm. 

The  plasma  density  profile  n{r)  and  the  thermal  energy  T are  obtained  by  dumping 
the  plasma  axially  and  measuring  the  total  charge  passing  through  a hole  in  a scanning 
collimator  plate.  Both  measurements  require  shot-to-shot  reproducibility  of  the  injected 
plasma,  and  we  typically  obtain  variability  bn/n  < 1%.  The  EV  plasmas  expand  radially 
towards  the  wall  with  a characteristic  “mobility”  time  of  t,„  100  sec,  so  the  electrons 

are  repetitively  injected  and  dumped.  On  IV,  a “rotating  wall”  perturbation  at  / ~ 
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FIGURE  2.  a)  Spectrum  of  me  = 0,  m^  = 1,  2,.. .5  Trivelpiece-Gould  modes  for  3 drive  amplitudes 
including  no  drive,  i.e.  thermally  excited. 


0.5  MHz  is  used  to  obtain  steady-state  confinement  of  the  electron  column  [6]. 

The  parallel  temperature  7j|  of  the  electrons  can  be  measured  by  slowly  lowering  the 
confinement  voltage  and  measuring  the  escaping  charge  [7,  8].  On  EV,  the  electrons 
equilibrate  to  T < 1 eV  soon  after  injection,  whereas  the  electrons  in  IV  cool  to  T w 
0.05  eV  due  to  cyclotron  radiation.  To  control  the  temperature,  we  apply  auxilliary 
“wiggle”  heating  by  modulating  one  electrode  voltage  at  a frequency  / = 0.8  — 1 .0  MHz. 
The  amplitude  and  duration  of  the  heating  is  adjusted  to  obtain  temperatures  up  to 
T < 5 eV;  higher  temperatures  cause  background  gas  ionization  and  a slow  increase 
in  N. 

We  perform  plasma  wave  transmission  experiments  by  applying  an  RF  voltage  of 
amplitude  V,,,  to  a (360°)  wall  electrode  at  one  end  of  the  column  (Fig.  la).  induces 
density  perturbations  5n  in  the  column  which  in  turn  induce  the  measured  voltages 
V,„  on  a (360°)  receiving  electrode  with  finite  load  impedance.  The  load  impedance  is 
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Rl  = 750  a (IV)  or  50  Q (EV)  in  parallel  with  Cl  ^ 440  pF  (IV)  and  Cl  = 165  pF  (EV). 

Figures  2a  and  2b  show  the  spectrum  of  standing  niQ  — 0 Trivelpiece-Gould  modes 
[9]  excited  with  excitations  of  V„,  = -80  dbm  (22  /tV)  and  -100  dBm  (2.2  jxV)  at 
frequencies  / = 0.01  - 10  MHz.  These  standing  waves  are  plasma  oscillations  with 
wavenumber  k-  = nm~/Lp.  Here  cOy,  = (2n)  28  MHz  («/10^  cm'“-^)'/“  is  the  plasma 
frequency,  modified  by  the  finite  trap  radius,  at  frequencies 


f ~ —r„in, — 
■’  In  ' -^L,, 


1 

1/2 

2 A../^p 

(I) 


The  axial  mode  number  is  = 1,2. ..5,  and  thermal  corrections  depend  on  the  ratio  of 
V = (kT to  the  wave  phase  velocity  v^.  In  this  paper,  we  focus  only  on  the  lowest 
radial  mode  nir  = 1. 

The  peak  amplitudes  are  proportional  to  V„.,  as  expected  for  a linear  .system;  the 
lesser  sensitivity  for  = 2 is  due  to  the  placement  Zc  and  length  (Lc  =11.7  cm) 
of  the  detection  cylinder.  The  peak  amplitudes  for  the  continuous  sinusoidal  modes 
are  independent  of  the  bandwidth  (BW  = 3 kHz)  of  the  spectrum  analyzer;  whereas 
the  spectral  amplitude  of  the  inter-mode  noise  decreases  as  (BW)"*/^  as  expected.  At 
V„,  = —100  dBm,  the  mode  fluctuations  have  amplitude  8n//i  10  The  peak  labelled 

RW  is  the  rotating  wall  drive;  the  measurements  presented  here  have  also  been  obtained 
with  the  drive  off. 

Small  peaks  representing  thermally  excited  modes  are  still  visible  in  Fig.  2c  when 
the  transmitter  electrode  is  grounded  (V,,,  = 0).  These  peaks  have  amplitudes  of 
-124  dBm/ \/3kHz,  representing  voltage  fluctuations  on  the  electrode  with  spectral 
intensity  F,,,  se  2.6nV'/v^. 

Figure  3 shows  received  spectra  of  the  thermally  excited  w-  = 1 mode  for  4 similar 
plasmas  of  different  temperature.  The  mode  frequency  fo  increases  with  temperature,  as 
expected  from  Eq.  (1).  The  width  of  the  spectral  peak  represents  mode  damping,  and  this 
width  increases  substantially  as  Landau  damping  becomes  significant  for  T > 0.5  eV. 

We  will  see  below  that  the  integral  / ^df  corrected  for  the  load  impedance  will 
correspond  to  the  thermal  energy  ^kT. 

Figure  lb  shows  a circuit  modeling  the  reception  of  thermal  noise  from  the  plasma. 
The  thermally  excited  fluctuating  plasma  voltage  V),  flows  through  a plasma  impedance 
Zp,  then  through  a load  impedance  Zl  with  its  inevitable  noise  V/.  Near  a mode  at 
frequency  o)o  = 2nfo,  the  plasma  admittance  Z''  is  given  by  a simple  pole,  as 


/■(to-o)())-fY,/ 


(2) 


As  will  be  seen  below,  the  plasma  resistance  Rp  depends  on  and  on  the  geometric 
coupling  between  the  plasma  and  the  receiving  electrode.  The  measured  load  resistance 
Rl  and  capacitance  Cl  give 


Z^  ' = R^  ^ -f  i(oCL 


(3) 


256 


Frequency  [MHz] 


FIGURE  3.  Spectra  of  the  thermally  excited  me  = 0,  = 1,  mode  for  different  plasma  temperature; 

the  solid  lines  are  fits  to  Eq.  (4). 


which  is  essentially  constant  over  the  mode  resonance;  this  simple  Zi  differs  from  the 
resonant-circuit  loads  commonly  used  in  harmonic  traps  with  a small  number  of  particles 
[3]. 

Nyquist’s  theorem  says  that  the  spectral  density  of  the  square  of  the  noise  voltage  is 
proportional  to  the  kT  times  real  part  of  the  impedance,  for  both  the  plasma  and  the 
load  noise  sources  [4],  A voltage-divider  fraction  Zj,/(Zp  + Z/,)  of  the  plasma  voltage 
Vp  will  be  measured  on  the  electrode  as  y„„  together  with  an  analogous  fraction  of  the 
(uncorrelated)  load  noise,  giving 


(4) 


where  = Re{Z}  and  Z*""  = Im{Z}.  Using  Eqs.  (2)  and  (3),  Eq.  (4)  can  be  explicitly 
written  as 


=4fcrz,z|^®  + 4 
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For  relatively  high  plasma  temperature,  one  can  assume  P = 0,  and  Eq.  (5)  describes 
a simple  Lorentzian  resonance  superimposed  on  a uniform  noise  background.  The  spec- 
trum is  then  characterized  by  4 parameters:  3 characterize  the  resonance,  as  frequency 
(Oo/2Jt,  width  Yiot.  and  amplitude  kT /Rp\  and  the  4th  parameter  of  the  load  tempera- 
ture Tl  sets  the  background  level.  Best-fit  values  of  kT /Rp  obtained  for  the  4 spectra  of 
Fig.  3,  together  with  the  separately  measured  plasma  temperature  T,  then  imply  values 
of  = 81,  160,  850  and  7800  Q. 

In  contrast,  both  T and  Rp  are  uniquely  determined  by  the  emission  spectrum  alone 
for  regimes  where  the  load-generated  noise  “filtered”  by  the  plasma  resonance  becomes 
significant,  as  in  the  T = 0.06  eV  case  of  Fig.  3.  Here,  the  spectrum  is  not  a simple 
Lorentzian,  and  a fifth  parameter,  Rp/Zf^,  is  required  for  a complete  fit  to  Eq.  (5).  The 
plasma  impedance  is  effectively  filtering  the  noise  generated  by  creating  a “dip”  in 
the  trace  at  coq-  The  imaginary  part  of  Z/.  creates  a small  shift  coq  - cOq  of  the  plasma 
peak  relative  to  the  noise  “dip.”  For  the  narrow  peak  of  T = 0.06  eV,  the  frequency 
shift  is  comparable  to  the  width  of  the  plasma  mode.  Under  the.se  conditions,  the  noi.se 
generated  by  is  effectively  performing  the  reflection  measurement  described  below. 
A 5-parameter  fit  to  this  spectrum  gives  T = 0.03  eV  and  Rp  = S\  Q,  consistent  with  the 
4-parameter  fit  and  separately  measured  temperature. 

For  most  of  our  data,  we  have  ncar-Lorentzian  V,„  spectra;  and  to  use  them  as  a tem- 
perature diagnostic,  Rp  must  be  obtained  independently.  This  plasma  resistance  Rp  (and 
full  impedance  Zp)  can  be  obtained  either  from  a reflection/absorption  measurement,  as 
shown  in  Fig.  Ic;  or  from  a first-principles  kinetic  analysis  of  the  plasma  fluctuations  [1] 
including  the  trap/antenna  geometry,  as  discussed  below. 

The  direct  reflection/absorption  measurement  of  Zp  uses  a directional  coupler  and 
lock-in  detector  to  determine  the  reflection  coefficient  r(f)  for  a weak  wave  at  frequency 
/ incident  on  the  receiving  antenna  and  plasma.  This  (complex)  reflection  coefficient 
is  defined  as  the  voltage  fraction  (and  pha.se)  which  is  reflected  by  the  plasma-loaded 
antenna  compared  to  that  reflected  by  an  open  circuit  without  antenna  or  plasma,  i.e. 


f = V'fcfl(plasma)/V'refl(open). 


(6) 
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FIGURE  4.  Measured  magnitude  and  phase  of  the  reflection  coefficient  f. 


This  reflection  coefficient  depends  on  the  impedance  Ztot  connected  to  the  directional 
coupler  compared  to  the  impedance  Zq  = 50  Q of  the  coupler  itself,  as  [4] 


Ztot  — Zq 


(1+r) 

(1-r)- 


(7) 


Here,  ZiJ  = Zp  ' + tcoC/,  is  the  total  impedance  of  the  plasma-loaded  antenna,  given 
by  the  plasma  impedance  Zp  in  parallel  with  the  capacitance  Cl  of  the  electrode  and 
connecting  cable. 

Figure  4 shows  the  measured  amplitude  and  phase  of  the  reflected  wave  as  the  fre- 
quency is  scanned  across  the  = 1 mode.  (The  data  is  from  EV,  where  the  frequencies 
are  about  2x  higher  than  on  IV.)  Away  from  the  mode,  the  entire  signal  is  reflected.  On 
resonance,  about  60%  of  the  incident  wave  is  absorbed  by  the  plasma,  and  40%  is  re- 
flected. A fit  to  the  data  gives  the  parameters  of  Z|ot  as  /o  = 4.37  MHz,  y/co  = 2.4  x 10“^, 
Rp  = 153  £2,  and  Cl  = 155  pF.  In  essence,  the  depth  of  the  absorption  dip  indicates  how 
close  Rp  is  to  the  50  Q of  the  directional  coupler. 

Alternately,  a kinetic  analysis  reproduces  the  Nyquist  theorem  and  gives  Rp  in  terms 
of  geometry  and  the  mode  damping  Yp-  Integrating  in  frequency  over  the  nie  = 0,  = 1, 

= 1 mode  gives  a mean-square  radial  electric  field  on  the  wall  of 


8E^  = kr 


8 


Rp  \ 

Rp  + Zl^) 


1 

4TtEo’ 


(8) 


with  J(x)  = \jI{x)  + j\{x)]l{dDldx)'^,  D{x)  = xJi{x)\n{rw/rp)  - Jo{x),  and 
X = (cOp/cOQ  — \y/^nrpjLp.  Here,  Jq  and  J\  are  Bessel  functions  of  the  first  kind, 
and  for  simplicity  we  have  kept  only  lowest-order  terms  in  rw/Lp.  Note  that  8Er  rep- 
resents the  radial  component  of  the  electric  field  at  the  wall,  not  the  total  wave  electric 
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Tdump 


FIGURE  5.  Plasma  temperature  measured  by  emission/reflection  technique,  compared  to  the  .standard 
dump  temperature  measurement.  The  triangles  are  from  the  IV  apparatus  and  the  diamonds  are  from  the 
EV  apparatus. 


field.  This  5£r  induces  the  measured  fluctuating  charge  5Q  on  a cylindrical  electrode  of 
length  Lc  located  at  given  by 


5Q^  = 2kT 


( 


ffp  ] 


g (4jieo) 


(9) 


with 


g = L^ 


sin 


^7tfe.  + Z.c)  j 


l2 


From  the  circuit  perspective  of  Fig.  lb,  0)q^Q^  represents  the  mode-integrated 
mean-square  current,  so  (0562“  = / df{Vp/df){Zp  + Zi)~^  —2kT  yi,/{Rp  + Zf'^)  using 
Eqs.  (2)  and  (3).  Thus,  we  obtain 


Ytot  = Yp  + Yz.  = ^ ) “0  g ) • (10) 


That  is,  the  kinetic  analysis  gives  mode-integrated  fluctuating  currents  proportional  to 
kT,  and  identifies  the  total  mode  damping  as  proportional  to  the  sum  of  the  resistances 
Rp  and  Zf^  in  the  plasma  and  the  load,  scaled  by  the  geometric  factors  g and  Ofo. 

For  typical  EV  parameters  of  Pp  — 1.66  cm,  r^^,  = 3.8  cm,  Lp  — 23.5  cm,  = 0, 
Lc  = 5.8  cm,  n = 8.45  x 10®  cm“^,  and  to^,/o)Q  = 41.3,  we  obtains  = 1.41,  D = 0.072, 
g = 1 cm/2.02  and  Y.ot/wo  - {Rp  + Zl^)/1.2  x lO^Q. 

It  is  important  to  note  that  contributes  to  the  mode  damping  on  an  equal  footing 
with  Rp.  For  example,  an  external  load  resistance  of  Zf^  — 12Q,  induces  a baseline 
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mode  damping  of  Yl/wq  = 10“^.  Thus,  even  small  receiver  resistance  can  contribute 
significantly  to  mode  damping. 

Figure  5 displays  the  inferred  plasma  temperature  Ter,  obtained  from  the  emission 
spectra  and  reflection/absorption  measurements,  versus  the  plasma  temperature  Taump 
measured  by  dumping  the  plasma.  The  plasma  temperature  was  controlled  by  adding 
“wiggle  heating”  at  one  end  of  the  column,  and  data  was  taken  for  plasmas  with  a range 
of  “geometric”  parameters  (n,  rp,  Lp)  on  both  EV  (diamond  symbols)  and  IV  (triangle 
symbols).  The  values  of  Ter  were  obtained  from  4-parameter  fits  to  the  emission  spectra, 
together  with  a separately  obtained  value  of  Rp. 

Conceptually,  determining  Rp  amounts  to  determining  the  plasma/antenna  geometry 
factor  ^ of  Eq.  (10).  Moreover,  Q does  not  depend  significantly  on  plasma  temperature, 
so  one  determination  of  ^ suffices  for  all  plasmas  with  the  same  geometric  parameters. 

On  EV,  we  measured  Rp  from  the  reflection/absorption  technique  over  a plasma 
temperature  range  of  0.5  <T  < 5eV.  This  gave  30  < Rp  < 2000£2,  implying  “internal” 
plasma  dampings  of  5 x 10“'*  < Yp/®o  < 0.02.  Over  this  temperature  range,  we  verified 
that  the  observed  Ytot  is  proportional  to  the  total  resistance,  with  a best-fit  value  giving 
^ = 1 cm/2.26.  This  is  in  close  agreement  with  the  algebraic  formula  for  Q. 

On  IV,  the  value  of  Rp  was  determined  from  a single  5-parameter  fit  to  the  non- 
Lorentzian  “T  = 0.06”  spectrum  of  Fig.  3,  giving  Q = 1 cm/5.3.  A direct  calculation 
of  Q using  the  “measured”  rp  = 0.2  cm  gives  Q = 1 cm/3.45;  but  this  discrepancy  is 
within  the  uncertainty  in  the  determination  of  small  plasma  radius  rp  on  IV. 

It  should  be  emphasized  that  steady-state  plasmas  often  exhibit  spectral  peaks  which 
are  10  to  100  times  larger  than  thermal,  because  the  mode  is  being  externally  driven. 
For  example,  noise  on  the  ostensibly  steady  confinement  voltages  Vc  or  ambient  RF 
signals  may  stimulate  particular  plasma  modes  without  proportionately  increasing  the 
plasma  temperature.  Similarly,  particular  modes  in  a warm  plasma  can  be  damped  (i.e. 
cooled)  with  an  external  feedback  circuit,  decreasing  the  peak  mode  amplitude  by  25  x , 
without  substantially  cooling  the  plasma.  That  is,  the  plasma  modes  are  relatively  weakly 
coupled  to  the  theorist’s  thermal  bath  [10].  Since  each  mode  has  only  kT  of  energy 
compared  to  the  10^  kT  in  the  entire  plasma,  it  is  rather  easy  for  particular  experimental 
manipulations  to  overwhelm  the  mode  energy;  so  this  thermometer  is  rather  fragile. 
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